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ABSTRACT 

Aims. We present results of three-dimensional nonlinear MHD simulations of a large-scale magnetic field and its evolution inside 
a barred galaxy with the back reaction of the magnetic field on the gas. The model does not consider the dynamo process. To 
compare our modeling results with observations, we construct maps of the high-frequency (Faraday-rotation-free) polarized radio 
emission on the basis of simulated magnetic fields. The model accounts for the effects of projection and the limited resolution of 
real observations. 

Methods. We performed 3D MHD numerical simulations of barred galaxies and polarization maps. 

Results. The main result is that the modeled magnetic field configurations resemble maps of the polarized intensity observed in 
barred galaxies. They exhibit polarization vectors along the bar and arms forming coherent structures similar to the observed ones. 
In the paper, we also explain the previously unsolved issue of discrepancy between the velocity and magnetic field configurations 
in this type of galaxies. The dynamical influence of the bar causes gas to form spiral waves that travel outwards. Each gaseous 
spiral arm is accompanied by a magnetic counterpart, which separates and survives in the inter-arm region. Because of a strong 
compression, shear of non-axisymmetric bar flows and differential rotation, the total energy of modeled magnetic field grows 
constantly, while the azimuthal flux grows slightly until 0.05 Gyr and then saturates. 
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1. Introduction 

A non-axisymmetric gravitational potential of the bar in barred 
galaxies influences strongly the dynamics of stars and the inter- 
stellar medium (ISM), including its magnetic field . The radio 
' polar ization observations of barred galaxies (e.g., iBeck et aP 
, 1200 5*) show that the magnetic field topology is far more compli- 
' cated than the gaseous component would indicate. No physical 
process responsible for such a complexity has been found so 
far . The radio maps of tvyenty b arred galaxies were pres ented 
■ in (iBeck et al.ll999ll2002ll2005l) and (iHarnett et al.l2004 . The 
' main observed features of t wo of the m ost spectacular galaxies 
NGC 1097 and NGC 1365 (iBeck et al 2005.) can be summa- 
rized as follows: 

I 

- the polarized emission (PI) forms ridges coinciding with 
dust lanes along leading edges of the bar in the central re- 
gion of galaxy; 

- the pitch angles of B-vectors vary rapidly in the bar region 
when located upstream from dust lanes, producing a depo- 
larization valley in the telescope beam; 

- the ridges of vanishing polarized intensity are observed in 
the shear shock areas; 

- in the outer disk, the polarization vectors form a spiral pat- 
tern with the maximum of emission along the spiral arms 
but in the inter-arm regions. 

The issue of magnetic field e volution in barred galaxies 
was a nalyzed in several papers (see Otmianows ka-Mazur et al.l 
120021 ll997HMosset alii 19981 [1999. 2001, 2007h . In Moss et 



al. (2007), three-dimensional (3D) calculations were applied to 
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study the process of magnetic field dynamo and gas evolution 
under the influence of a bar potential. Moreover, they prepared 
the global synthetic polarization maps to compare them di- 
rectly with the NGC 1365 radio observations. Even though all 
of the input parameters have been studied varying each inturn 
failed to reproduce a configuration of polarization vectors sim- 
ilar to that observed. The idea of magnetic arms betw e en the 
gaseous counterparts was also studied by Rohde etakl (Il998l) 
using 3D MHD dynamo simulations. The authors appUed a 
non-uniform distribution of magnetic diffusivity with its higher 
values in spiral arms and lower in the inter-arm region. This ar- 
rangement was related to the presence of higher turbulent mo- 
tions in the spiral arms due to the presence of star-forming ac- 
tivity there. The assumption of higher magnetic diffusion in 
spiral arms, no matter how reasonable, has never been con- 
firmed by observations. The authors found that it was possi- 
ble to obtain magnetic arms in the inter-arm region in the case 
of resonance of the density pattern speed and a magnetic drift 
velocity, thus not for all galaxies. Our nonlinear model of mag- 
netic field evolution in barred galaxies, which takes into ac- 
count the back reaction of magnetic field on the gas, offers an 
opportunity to understand the relationship between the large- 
scale magnetic and gas velocity fields. 

Our previous 3D simulations (lOtmianowska-Mazur et alJ 

120021) concerning the large-scale magnetic field evolution in 
barred galaxies applied the evolution in velocity fields pre- 
dicted by the N-body sticky-particle simulations at that time. 
It was impossible to introduce either the back-reaction of mag- 
netic field on the gas component or the physical mechanism 
of gaseous spiral waves on the magnetic field structure. In the 
present paper, we study the nonlinear evolution of the gas and 
magnetic field in a barred galaxy. Our 3D calculations are car- 
ried out with the use of a MHD code written by Kowal (2008). 
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We do not include the dynamo process in our MHD equations. 
The gravitational potential originating in stars in the form of the 
bar, halo, and bulge is modeled analytically. The magnetic field 
of galaxies is usually studied using polarized radio-continuum 
observations, yielding the magnetic field configuration inte- 
grated along the line of sight projected onto the sky plane and 
convolved to a certain radio telescope beam. Since our models 
contain full three-dimensional magnetic field structures, we an- 
alyze them by simulating the polarization maps and comparing 
these maps directly with observations. 

2. Methods 

2. 1 . The initial conditions and input parameters 

We investigate the evol ution of a barred galaxy s olving the re- 
sistive MHD equations (iLandau & Lifshit3l998l) . We apply an 
isothermal equation of state pgas = pc^, where c., - 5 km s ' is 
a constant isothermal speed of sound. In our model, we assume 
that the initial gas density distribution is of the form: 

p=poexpJ-Mj2(r), (1) 

where po = 1 H cm"^ is the initial gas density of the neutral 
hydrogen in the galactic mid-plane, h = 0.3 kpc is the gas scale- 
height, and function Q(r) is the truncation factor that reduces 
the density to small values at r > R,„ax, Rmax - 9 kpc being the 
galaxy radius. We assume the azimuthal initial magnetic field 
configuration (B~ - Q, Br - 0, B^{z, r)). Its distribution de- 
pends strictly on the gas density distribution via the condition 
B = where B - 10.0 at the beginning (determining the 
initial value of magnetic field to be Z?o = 1 .0/vG in the galactic 
mid-plane). We use the turbulent dijfusion coefficient for the 
interstellar gas ?; = 3 ■ 10^^ cm^ s"' ( lLeschHl993l) . 

The gravitational potential of stars is introduced analyti- 
cally. Our galaxy is initially composed of three components: 
the large and massive halo, the central bulge and the rotating 
disk of stars. The circular speed of the gas is: 



+ GMh(r^ + alr^l^ + GMh{r^ + alY^'^ + 
c]{\+l3)r dQ 

where Md - 10'" M© is the mass of the disk of stars, Ud - 0.2 
kpc is its characteristic scale-length, and both Mh - 10"^ Mq, 
and Gb - 2 kpc refer to the bulge and M/, - 10" Mq, and 
fl/i - 20 kpc refer to the halo. The bar component grows grad- 
ually in time (from theg = O.lGyr to te„d - OAGyr), chang- 
ing the gravitational potential of the galaxy. The bar poten- 
tial i s described by the second-order Ferrers ellipsoid dFerrersI 
11877 ) with semi-axes a - 4kpc, b - 2 kpc, and c - 2 kpc. 
To conserve the total mass of the galaxy, we reduce the bulge 
mass, such that M^arit) + Mi,{t) - M, where M is constant 
during the calculations. The bar angular velocity is set to 
be 30kms ' kpc This determines the value of the corotation 
radius to be Rcor = 4.46 kpc. 

We use a higher-order shock-capturing Godunov-type 
scheme based on the essentially non-oscillatory (ENO) spa- 
cial reconstruction, a Harten-Lax-van-Leer (HLL) approxi- 
mate Ri^nwnn;Soh[e£_and^ time integra- 
tion (iDel Zanna. Bucciantini & Londrilldl2003L e.g.) to solve 



isothermal non-ideal MHD equations. We evolve the induc- 
tion equation using t he constrained transport (CT) method 
(lEvans. Hawleyll988h to ensure that magnetic field divergence 
vanishes everywhere at all times. 

The computational domain extends from - 10 kpc to 10 kpc 
in the x and y directions, and from -2.5 kpc to 2.5 kpc in the z 
direction with the grid size rix = riy = 512, - 129. 

2.2. Stokes parameters 

To compare our results with observed properties of barred 
galaxies, we constructed synthetic radio-polarization maps 
(face on) at selected time steps. We calculated the maps by in- 
tegrating Stokes parameters of the synchrotron emissivities I, 
Q, and U along the line of sight and finally convolving them 
with an assumed Gaussian beam. As a result, we generated po- 
larization maps of the intensity and the angle of polarization 
B-vectors. We assumed that the distribution of relativistic elec- 
trons change radially and vertically as Gaussian functions. The 
truncation radius and scale height were respectively equal to 
9 kpc and 1 kpc. The adopted intrinsic polarization degree of 
synchrotron emission is 70% and the energy spectral index of 
the relativistic electrons y is 2.8. 



3. Results 

3. 1 . Velocity structures and magnetic field amplification 

Figure [1] shows the gas density distribution (grey-scale plot) 
and velocity field (in a reference frame corotating with the bar) 
in the galactic mid-plane at time t = 0.64 Gyr. The bar's ma- 
jor axis is inclined by 20° to the x-axis. The high density re- 
gions appearing from the innermost region of the galaxy in- 
dicate the leading edge of the bar Density enhancements are 
visible at the ends of the bar and in the spiral arms generated 
in the galactic disk. The bar strongly disturbs the gas velocity 
field. Rapid variations in the velocity vector direction can be 
seen in the innermost region of the bar The gas flows out of 
the center just before it passes the shock, and then turns back 
toward the center along the leading edge of the bar. Gas mo- 
tions generate shocks that strongly influence the magnetic field 
structure. The bar is gradually introduced to the model from 
tbeg = 0.1 Gyr to te„d - 0.4 Gyr. During this time, as we can 
see in (Fig. |2]l, the mean value of the B^ flux in the galactic 
mid-plane increases slightly and then saturates. No significant 
enhancement in the B^ flux can be seen. The total magnetic 
field energy (B^) increases as the galaxy evolves, because of 
an increase in the magnetic field complexity, which is in turn 
the result of a complex velocity field global structure. In our 
model, the most significant part of the growth in the magnetic 
field energy originates in the combined effects of shearing mo- 
tions and compression present in the innermost bar region. 

3.2. The polarization maps 

Figure [3] shows the distributions of polarization angle (vectors) 
and polarized intensity (contours) superimposed on the col- 
umn density (grey-scale) at time t=0.64Gyr. To illustrate the 
evolution of the magnetic field configuration in our model, we 
present results for three crucial time steps: 0.3 Gyr (Fig. 3 left), 
0.46 Gyr (Fig. 3 in the middle), and 0.64 Gyr (Fig. 3 right). As 
mentioned above, in the beginning the modeled barred galaxy 
is in magneto-hydrostatic equilibrium. After almost one and a 
half rotation periods (fi - 0.30 Gyr), the density distribution 
diverges significantly from its initial form. In the inner part of 
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Fig. 1. Velocity field in the rest frame of the bar (vectors) and 
the gas density distribution (grey-scale, dark means high values 
of density) along the z = plane at time t - 0.64 Gyr. 



the disk, where the bar is present, we find the highest gaseous 
density and magnetic field intensity regions. 

The compression of both the gas and magnetic field is 
also visible in the outer part of the disk in the form of spiral 
arms. All of these gaseous structures are created by the dy- 
namical influence of the non-axisymmetric gravitational po- 
tential of the bar. This potential also changes the initial az- 
imuthal configuration of the magnetic field. First, the magnetic 
field follows the gas distribution, as can easily be seen for time 
step t\ — 0.30 Gyr in Fig. [3] (left), where the magnetic field 
strength maxima are aligned with the gaseous ones. However, 
at later time-steps, the magnetic arms begin to detach them- 
selves from the gaseous spirals and drift into the inter-arm re- 
gions preserving similar to the gaseous arms pitch angles (see 
f2 = 0.46 Gyr Fig. |3] middle). At f3 = 0.64 Gyr (Fig. |fl right), 
we can observe that our magnetic arms are distributed between 
density wave structures. The magnetic arms do not corotate 
with gaseous spiral structure, but they follow the general gas 
motion in the disk, which has a slightly lower angular velocity. 
The last map (Fig. [3] right) also shows the depolarization re- 
gions present in the bar. They are the effect of so-called beam 
depolarization, where magnetic field lines both along the bar 
and in the inter-arm regions converge with large angles within 
s mall, one-beam- size areas. Similar results were also obtained 
by lOtmianowska-Mazur et al ] (l2002h . 

3.3. Magnetic pitch angles 

To analyze and compare the pitch angles of the magnetic and 
gaseous arms, we project our results onto the polar coordi- 
nate space of azimuthal angle in the disk and ln(r), where r 
is the galactocentric distance. In this case, the logarithmic spi- 
ral is represented by a straight line inclined by its pitch angle. 
The grey plot of the gas density (integrated along the line of 
sight) and superimposed contours of polarized intensity and B- 
vectors at a time 0.64 Gyr are shown in Fig. |4] The values of 
the pitch angle of the magnetic arms are similar to those of the 
gaseous ones, i.e., they are of order of 5°. We also estimated the 
mean pitch angle (averaged over the azimuthal angle and radius 
in the plane of galaxy) to be 3°. Since we do not apply dy- 
namo effects, the aforementioned pitch angle values are rather 
low, slightly smaller than usually observed in spiral galaxies 
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Fig. 2. Time dependence of the total magnetic field energy 
(solid line) and the mean B^ flux in the galactic mid-plane 
(dashed). The third plot (dotted) presents the total magnetic 
field energy without taking into account the innermost bar re- 
gion. All values are normalized to their initial values. 



(see lBecklll993h . However, between the radii r - 4.1 kpc and 
r = 4.5 kpc, the mean pitch angles are much higher and reach 
24°. 

In the inner part of the disk (between the bar and spiral 
arms), the B-vectors of the polarized intensity change their 
pitch angles rapidly (Fig H] see <p around 100 and 300, and 
In r between 2.2 and 3.2). This is caused by the strong flow 
of the gas between the end of the bar and the magnetic spiral 
arms, which produces valleys of depolarization within a tele- 
scope beam. T his bears some res emblance to the observations 
of NGC 1356 (lBecketal.ll2005h and NGC 3627 (ISoida et alJ 
1200 Ih , even though we applied different rotation curves. 



Fig. 4. The distribution of polarized intensity and B- vector ori- 
entation superimposed on the gas density at /i = 6.2 cm in the 
frame of azimuthal angle in the disk and the logarithm of dis- 
tance from the center. The map has been convolved to the res- 
olution 40". 



4. Discussion & Conclusions 

For the first time, we have simulated the evolution of large- 
scale galactic magnetic fields in baiTed galaxies using 3D MHD 
nonlinear numerical simulations without a dynamo process. 
The physical processes modeled and observed in our simula- 
tions can be summarized as follows. First, the dynamical influ- 
ence of the bar causes the gas to form spiral waves. Later on, 
both the shear and compression in the arms cause the magnetic 
arms to be aligned with the gaseous structures. Since the pat- 
tern of the gaseous spiral arms travels faster than the disk itself, 
especially beyond the corotation radius, the gaseous structures 
move further out, while the magnetic arms remain in the inter- 
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Fig. 3. Face-on polarization maps at A = 6.2 cm at selected times steps (0.3 Gyr left, 0.46 Gyr middle, 0.64 Gyr right): polarized 
intensity (contours), polarization angle (vectors) and column density (grey-scale). AU maps have been smoothed down to the 
resolution 40". The black color represents the regions with the highest density. 



arm regions. This results in the magnetic and spiral arms having Rohde, R., Elstner, D. & Rudiger, G., 1998, A&A, 329, 911 

similar pitch angles. Soida, M., Urbanik, M., Beck, R., Wielebinski, R., Balkowski, C, 2001, A&A 

We found that: 



1. Magnetic arms develop inside the gaseous ones, but grad- 
ually diverge away from the desnity waves. The magnetic 
arms are consequently shifted toward the interarm regions, 
in qualitative agreement with observations. 

2. The above effect is also responsible for prducing the de- 
polarized valleys in the bar region, which also agrees with 
observations. 

3. The magnetic pitch angles of the synthetic magnetic arms 
are similar to the pitch angles measured in spiral gaseous 
structures, which also agrees with observations. 

This process could also be responsible for the shifting of 
the magnetic arms to the interarm regions in spiral galaxies. 
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